UNIQUENESS OF TRANSONIC SHOCK SOLUTIONS IN A 
DUCT FOR STEADY POTENTIAL FLOW 
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Abstract. We study the uniqueness of solutions with a transonic shock 
in a duct in a class of transonic shock solutions, which are not necessar- 
ily small perturbations of the background solution, for steady potential 
flow. We prove that, for given uniform supersonic upstream flow in a 
straight duct, there exists a unique uniform pressure at the exit of the 
duct such that a transonic shock solution exists in the duct, which is 
unique modulo translation. For any other given uniform pressure at 
the exit, there exists no transonic shock solution in the duct. This 
is equivalent to establishing a uniqueness theorem for a free boundary 
problem of a partial differential equation of second order in a bounded 
or unbounded duct. The proof is based on the maximum/comparison 
principle and a judicious choice of special transonic shock solutions as a 
comparison solution. 



1. Introduction and Main Results 

In the recent years, there has been an increasing interest in the study 
of transonic shock solutions in ducts or nozzles for the steady potential 
flow equation or steady full Euler system for compressible fluids. The basic 
strategy is to construct first some transonic shock solutions and then study 
the stability of these solutions by perturbations of the boundary conditions; 
see [U 13 [3 [15j [HI [18] and the references cited therein. On the other 
hand, some basic properties of these special transonic shock solutions, such 
as the uniqueness in a large class of solutions, have not been fully understood. 
As a first step, in this paper, we study the uniqueness of solutions with a 
flat transonic shock in a straight duct in a class of transonic shock solutions, 
which are not necessarily small perturbations of the background solution, 
for steady potential flows. Some classical, related results on transonic flows 
may be found in (9j [13] and the references cited therein. 

Consider steady isentropic irrotational inviscid flows in a finite duct D := 
(—1,1) x C K 3 or a semi-infinitely long duct D' := ( — 1, oo) x VL C M 3 , 
where C I 2 is a bounded domain with C 3 boundary. 
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The governing equations of potential flows are the conservation of mass 
and the Bernoulli law (cf. |10j): 

V • (pVip) = 0, (1.1) 

±\V<p\ 2 + i(p) = b , (1.2) 

where (p is the velocity potential (i.e., V</? is the velocity), 6o is the Bernoulli 
constant determined by the incoming flow and/or boundary conditions, p is 
the density, and 

p p 

with c(p) being the sound speed and p(p) the pressure. For polytropic gas, 
by scaling, 

p(p) = ^-, c\p) = p<-\ i( p ) = ^—Zl, 7 >1. (1.3) 
7 7-1 

In particular, when 7 = 1 as the limiting case 7 — ► 1, 

i{p)=hxp. (1.4) 

Expressing p in terms of |V<^| 2 : 

p = p(\V V \ 2 ) = (l + (7 " l)(6o " \\^^\ 2 ))~ 1 for 7 > 1, 

or 

p = p(V^\ 2 ) = e -iiv^l 2 +b fo r7 = i ) 
equation (II. ID becomes 

V • (p(\Vip\ 2 )Vip) = 0. (1.5) 

Equation (|1.5p is a second order equation of mixed elliptic-hyperbolic type 
for ip in general; it is elliptic if and only if the flow is subsonic, i.e., |Vy| < c 

or equivalently, \V<p\ < c* := ^/ (l + (7 — l)&o) for 7 > 1 and c* = 1 for 
7 = 1. 

We first consider the case of a finite duct D. Let V = [—1, 1] x <9f2 be the 
lateral wall, and let Sj = {i} xfl, i= —1, 1, be respectively the entry and 
exit of D. That is, we assume 

(ill) d <p>0 on £< = {«} xft, i = -1,1. 

Note that dD = E_i U £1 U T. We are interested in the case that the flow is 
uniform and supersonic (i.e., |Vy| > c) on S_ 1; subsonic on Si with uniform 

pressure. More specifically, for a constant u~ G (c*, \/2(bo + ^1)) and a 



UNIQUENESS OF TRANSONIC SHOCK SOLUTIONS 



3 



constant c\ G (0, c*), we consider the following problem: 

CCS) in A (1.6) 

9? = — u~ , d Xl ip = u~ on (1-7) 

|V^| = ci on Ei, (1.8) 

V(^-n = on T, (1.9) 

where n is the outward unit normal on T. 

We remark that the formulation of this boundary problem is physically 
natural. Since the flow is supersonic near £_i, i.e., the equation is hyperbolic 
on £_i, there should be initial data like (|1.7p due to (-Hi) (Our choice 
of ip in (jl.7p makes the solution of the uniform supersonic upstream flow 
in D looks neatly; see Lemma 11.11 below). On the other hand, since the 
equation is elliptic on £i, only one boundary condition is necessary. We 
choose the Bernoulli-type condition (jl.8p since, from the physical point of 
view, assigning the pressure (i.e. density for isentropic flow) is of more 
interest (cf. [10j). which is just a boundary condition like (jl.8p due to 
(jl.2p . Condition (jl.9p is the natural impermeability condition, i.e., the slip 
boundary condition, on the lateral wall for inviscid flow. 

We are interested in the class of piecewise smooth solutions with a tran- 
sonic shock for problem (jl.6j) - (jl.9p . 

Definition 1.1. For a C 1 function x\ = f(x 2 ,x 3 ) defined on Q, let 

S = {(f(x 2 , x 3 ),x 2 ,x 3 ) G D : (x 2 , x 3 ) G £1}, 

D~ = {(x 1 ,x 2 ,x 3 ) G D : x\ < f(x 2 ,x 3 )}, 

D + = {(x 1 ,x 2 ,x 3 ) G D : x\ > f(x 2 ,x 3 )}. 

Then ip G C 0,1 (-D) D C 2 (D~ U D + ) is a transonic shock solution of (jl.6D 
(I1.9P if it is supersonic in D~ and subsonic in D + , satisfies equation (II. 5p in 
D~ U D + and the boundary conditions (|1.7j )- (|1.9j) pointwise, the Rankine- 
Hugoniot jump condition: 

p(\Vip + \ 2 )X7tp + -v = p(\Vip-\ 2 )V<p- -v onS, (1.10) 

and the physical entropy condition: 

KIV^+I 2 ) > p(\V^-\ 2 ) ^ \V(p + \ <\Vip~\ on S, (1.11) 

where v is the normal vector of S, and ip + (<p~) is the right (left) limit of (p 
along S. The surface S is also called a shock-front. 

Remark 1.1. Note that, across a transonic shock-front S", the potential 99 is 
continuous, while the velocity V99 is discontinuous. Since the shock-front 
is a free boundary that requires to be solved simultaneously with the flow 
behind it, we have to deal with a free boundary problem indeed. In the 
following, we also write (p^ = p\o ± with ip~ the supersonic flow and (p + the 
subsonic flow. 



We first state two direct facts. 
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Lemma 1.1. There exists a unique supersonic flow solution tp~ that satisfies 
(|1.6p ~ (|1.7j) and (11. 9ft in i/ie c/ass o/C 2 supersonic flow solutions in D. The 
unique solution is (p~ = u~x\. 

Proof. It is clearly that (p~ = u~x\ solves (jl.6|) - (jl.7p and (jl.9p . By standard 
energy estimates for hyperbolic equations, this solution is unique in the class 
of C 2 supersonic flow. □ 

Lemma 1.2. For each t € ( — 1, 1), the function with a flat transonic shock- 
front: 

, \ juTxi, -1 < XX < t, 

ipt{x 1 ,x 2 ,x 3 ) = < (1.12) 
lr(xi- t) + u t, t < xi, 

solves problem (jl.6)) (|1.9f) with c\ = u + , where u + £ (0, c*) is determined 
by u~ , bo, and 7 > 1. 

Proof. This is equivalent to solving u + from the following two algebraic 
equations deduced from (II. 2\ and (11,101) : 

p~ u - = p + u + , (1.13) 

iq^ + ^?=^ f ± + i (u - r . (1 , 4) 

7— 1 z 7 — 1 z 

The calculation similar to Proposition 3 in [16] indicates that there exists 
a unique solution u + < c* < One can then easily verify that ipt con- 
structed above is a transonic shock solution. □ 

We remark that this special solution has played a significant role in the 
study of transonic shocks in the recent years for the potential flow equation 
and the full Euler system (cf. d EJ EJ [HJ 13 Q21 IB] ) . It has been observed 
to be unstable if the pressure is given at the exit in general. Theorem 11.11 
below provides a simple and direct confirmation of this instability. 



Theorem 1.1. Under assumption (Hi), for given u E (c*, w 2(6o + ^hf))> 

problem (jl.6p ~ (|1.9p is solvable for transonic shock solutions in the sense of 
Definition 1.1 if and only if c\ = u + , with u + being a constant in (0, c*) 
determined by u~ , bo, and 7 > 1. In addition, the solution is unique modulo 
translation: it is exactly ipt for t € ( — 1, 1). 

This result especially implies that, if the pressure is posed at the exit, 
then the boundary value problem is ill-posed in most cases, and the spe- 
cial transonic shock solutions (p t that are widely studied are not physi- 
cally stable (cf. [U [HI E])- We remark that, unlike the previous works 
[H HJ El [71 HJ [TSl HTl JS] where the stability of the special solutions was 
studied under small perturbations of the upstream supersonic flow (jl.7j) or 
the shape of the wall T of the duct, our results do not require such small 
perturbations. Our proof is global and based on the maximum/comparison 
principle and a judicious choice of special transonic shock solutions as a 
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comparison solution. It reveals the basic uniqueness property of such spe- 
cial transonic shock solutions. 

Next, we focus on the uniqueness of transonic shock solutions in semi- 
infinite duct D' = (—1, oo) x fi, with the following assumption only for the 
case 7 = 1 that 

(H2) there exists ko > such that |Vy| is bounded in (ko,oo) x 17. 

We remark that this assumption is automatically satisfied for potential 
flow with 7 > 1, since the velocity should be less than the critical value 

Y^2(6o + ^rr) due to the Bernoulli law and the fact that the constant 60 has 

been fixed by the supersonic data at the entry. 

Let r = (—1, 00) x dd. We have the following result: 

Theorem 1.2. Consider problem (|1.6|) - (|1.7p and (jl.9p with D replaced by 
D 1 and (II. 8p replaced by (H2). Assume that, for the solution if, \D 2 (p\ is 
also bounded in (ko,oo) x Q. Then this problem is solvable for transonic 
shock solutions, and the solution is unique modulo translation. 

Remark 1.2. This result indicates that, for transonic shock solutions in semi- 
infinitely long ducts, there should be no additional asymptotic condition 
such as 

lim max \V<p(xx, x%, x 3 )\ = c%. (1-15) 

Otherwise, it is either overdetermined or superfluous. We just need the 
reasonable assumptions that the velocity and acceleration of the flow are 
bounded. This indicates that the apriori assumptions on the asymptotic 
behavior of transonic shock solutions in an infinitely long duct in [5j El [7] 
may not be necessary. 

Theorem 1.3. In Theorem II. 21 if we replace (H2) by the following stronger 
assumption: 

{H'2) I Vy| < c < c* in (ko, 00) x £1 with c a constant, 

that is, the far-away- flow field behind the shock-front is always subsonic, 
then the requirement that \D 2 f\ is bounded can be removed. 

In the rest of this paper, Sections 2-4, we establish Theorems I1.1H1.31 
respectively. 

2. Proof of Theorem 11.11 
We divide the proof into three steps. 

Step 1. Let cp be a transonic shock solution of problem (|1.6jl - (|1.9jl . and 
let 5 be the corresponding shock-front with equation x\ = f(x2,xs), and 
t = mixin{f(x 2 ,x 3 )}. Then 

ip = u~f(x 2 ,x 3 ) on S, (2.1) 
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since the potential <p is continuous across the shock- front. 

Step 2. We notice that (12, lj) and the Neumann condition (|1.9p imply that 
S is perpendicular to V. In fact, for P G T n 5", let the normal vector of 
r at P be n = (0,n2,n3), and denote the normal vector of S at P to be 
v = (1, —d X2 f, —d X3 f). We now show that re • v = 0. 

By (|2.ip . we obtain that, at P, 

d Xi ip + d Xl ipd x J = u~d x J for i = 2,3. (2.2) 
Therefore, by Ojl . 

3 3 

(it~ - d xi ip) y^riid Xi f = yinjdxjipi = 0. (2.3) 
i=2 i=2 

By (jl.lip . after passing S 1 , we have < u~. Therefore, Yli=2 n i^if = 
and n ■ v = 0. 

5fep 5. Now let (p a := ip T be the transonic shock solution constructed in 
Lemma ll.2( D + and D£ be the corresponding subsonic region of <p and <p a , 
and D* := n -D + = D + . Then ip = (p a — <p satisfies the linear, uniformly 
elliptic equation: 

Lip = } y aij(x)d XiXj ip + bi(x)d Xi ip 

:= p(|V93 a | 2 )A^ - p'(\Vip a \ 2 )d Xi ip a d Xj ip a d XiXj i) 

+&p{p(\Vp a \ 2 )-p(\Vp\ 2 )) 
-2 d XiX . p (p'(\V(p a \ 2 )d Xi (p a d Xj p a - p'{\Vip\ 2 )d Xi ipd Xj ip\ 

= in D*. (2.4) 

Since D is bounded, by our assumption ip € C 2 , \V 2 p\ is bounded in D + . 
The boundary conditions are 

V(p a + p)-Vij = {u+) 2 -c 2 on S 1; (2.5) 

VV>-n = onmF. (2.6) 

By they are both the oblique derivative conditions. The boundary 

condition on S* = {(x±, x%, X3) : x\ = max{r, f(x2, X3)} = S is 

V> = g(x 2 , x 3 ) := (u+ - u-){f{x 2 ,x 3 ) - r) < 0. (2.7) 

Note that there exists Y G O such that f(Y) = r, so 

iP(f(Y),Y)=g(Y)=0. (2.8) 

We now prove that ip = in D*. By (|2.8p . it suffices to show that ip is a 
constant. There are two cases. 

Case A. u + > c\. By the strong maximum principle, m = min^* ip can 
be achieved only on dD* unless ip is a constant. 
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By the Hopf lemma (cf. Lemma 3.4 in [12j). the minimum m of ip can be 
achieved only on S or T D Si, but not the lateral boundary T and the exit 
Si unless ip is a constant. 

(i) Suppose that m is achieved at a point P £ T n Si. By a locally even 
reflection with respect to V and noting that V is perpendicular to Si at P, 
P satisfies the interior sphere condition in the extended neighborhood, as 
well as d Xl ip > due to (|2.5p and (|2.6p . a contradiction to the Hopf lemma 
unless ip is a constant. 

(ii) Suppose that m is achieved on S. Then m < 0. 

(a) . Let m = for some lefi. Note that, by ([221), V/(X) = 0, so 
v = (1,0,0). By the Rankine-Hugoniot condition (jl.lOp and the Bernoulli 
law (jl.2|) . as in Lemma ll.2l we can solve that 

d x MX)=d xl ip a (X)=u+. (2.9) 

Hence, 

Vi>(f(X),X)-v(f(X),X)=0. (2.10) 
By the Hopf lemma, it is impossible unless ip is constant. 

(b) . Let m = g(X) for some X € dfl. Then it is still necessary to hold 
V/(X) = due to the orthogonality of S and T. We also need a locally 
reflection argument as in (i) to apply the Hopf lemma to infer that ip is a 
constant as in (a) by (|2.10j) . 

Case B. u + < c\. Similar to the analysis in Case A, now the maximum 
M of ip in D* can be achieved only on S unless tp is a constant. According 
to (HZD, M = and we may also obtain Vip(f(Y),Y) ■ u(f{Y),Y) = 0, a 
contradiction to the Hopf lemma unless ip is constant. 

Therefore, ip = 0. This implies that, for c\ ^ u + , there is no solution; 
for ci = u + , the solution (p T is unique (i.e., for any given r G (—1,1), 
there is only one transonic shock with its front passing a point in {r} x Q, 
which is exactly <p T ). Since (pt+ T (xi, £3) = <^ T {xi — t,X2,xs), then, for 
c\ = n + , the solution to problem (jl.6p - (jl.9p is unique modulo translation 
in the x\ -direction. 

This completes the proof. 

Remark 2.1. Note that, in the proof, we need the assumption only that the 
downstream flow on the shock- front is subsonic (see (|2.9p and (|2.3p above). 
We do not need to assume that the flow behind the transonic shock-front 
is always subsonic. In addition, the assumption c\ < c* is needed just to 
guarantee that the transonic shock-front is restricted in D. 

3. Proof of Theorem 11.21 
We divide the proof into three steps. 

Step 1. Let <p be a transonic shock solution of problem (ll.6p - (|1.9p with 
D replaced by D' and (II. 8p replaced by (^2), as well as \D 2 tp\ is bounded 
in its subsonic region. Denote its shock-front as S = {(/(X2, X3), X2, X3) : 
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(x2,xs) € Q}. Let r = min^ / and ip = (p T — ip, where (p T is the special 
transonic shock solution constructed in Lemma 11.21 The key point is to 
show either the maximum or the minimum of ip in D'^ := {(#1, X2, X3) £ 
D' : x\ > f(x2,xs)} is achieved on S. With this, then the rest of proof 
is the same as that for Theorem 11.11 To achieve this, we now show the 
following case, Case A, is impossible if ip is not a constant. 

Step 2. Case A: Both the maximum (might be 00) and the minimum 
(might be —00) are achieved as x\ — > 00. 

We deduce below that, for this case, there is a contradiction if ip is not a 
constant. 

Define D' L = {(xi,X2,xs) € D' : f(x2,xs) < xi < L}, for L large, such 
that S C D' L . Without loss of generality, we assume Mi = max^ ip > 
0,mi = min^ ip < (cf. (12. Tf) ) . since, for the case Mi = or mj, = 0, we 
can still apply the Hopf lemma to show ip = 0. 

Note that Mi is monotonically increasing, while mi is monotonically 
decreasing, as L increases. In this case (Case A), by assumption, for large 
L, both Mi and mi are achieved on T,i = {L} x unless ip is constant 
(which is what we want to prove). Since f2 is bounded and | Vip\ < |V<£>| +u + 
is bounded according to (#2), we conclude that both m = lim^oo mi < 
and M = lim/^oo Ml > are finite. Thus, ip is a bounded solution. 

Now choose a sequence {Lk}'j£ =1 that tends to infinity. On S fc := S_L fe , 
we suppose M k := M Lk = ip(L k ,X k ) > 0, m k := m Lk = ip(L k ,Y k ) < 
for X k ,Y k S $7. Since (p is continuous, there exists Z k € 0, such that 
ip(L k ,Z k ) = 0. 

The following arguments are adopted from [2]. Consider B k := (L k — 
2,L k + 2) x Q in M 3 , by suitable translation, each B k (k € Z) may be 
transformed onto B := (—2,2) x Q. 

Let ip k (Y) = ip((L k ,0) + Y), k = &q, ko + 1, • • • , which is defined on B and 
satisfies the linear elliptic equation: 



where a>*MY) = a i:j ((L k ,0) + Y) and bf\Y) = bi({L k ,0) + Y). Obviously, 



the ellipticity constants of these equations are the same, and the coefficients 
are also uniformly bounded. 



Now let K = (—1, 1) x Q be a relatively open set of B U ((—2, 2) x V) and 
K C BU ((-2,2) x T). 

Consider v k {Y) = M — ip k (Y), which is positive by definition of M. By 
(|3.ip . we have 



3 



3 




(3.1) 



3 



3 




(3.2) 
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Applying the boundary Harnack inequality for the oblique derivative prob- 
lems (cf. pQ, Theorem 2.1, or [14]) to Vk, we have 

C(M - MYi)) < Csupv k < inivk < (M - ^ k {Y 2 )) (3.3) 

K K 

for any Y\, I2 £ and the positive constant C is independent of k. Taking 
Y\ = (0, Zk) and Y2 = (0, X k ), and letting fc — > 00, we obtain 

CM < 0, (3.4) 

which is a contradiction to our assumption that M > 0. 
Therefore, Case A is impossible if ip is not constant. 

5fep 3. Suppose that ip is not a constant. By the strong maximum prin- 
ciple applied in the domain D'^, at least one of the maximum and minimum 
should be achieved on S. Applying the Hopf lemma on S as in Section 2, 
we may also infer ip = 0. This also contradicts our assumption that ip is not 
a constant. 

Therefore, we conclude that tp is constant. By our choice of r, it should 
be in D'^ . This completes the proof of Theorem 11.21 

4. Proof of Theorem 1.3 

We now show that, for a C 2 transonic shock solution ip to equations (|l.ljl - 
(jl.2p . if (H' 2 ) holds, then \D 2 ip\ is bounded in the unbounded subsonic region. 
The proof can be achieved by standard elliptic arguments as sketched below. 

Step 1. Decomposition of unbounded domain. For k S Z (positive inte- 
gers), let 

D k = (k-i,k + 2)xn, D' k = (k- ~ k + ^) x n. 

Clearly, we have dist(-D&, D' k ) = ^. 

Step 2. Uniformly boundary Holder estimate of gradient and second order 
derivatives of ip. 

For s = 1, 2, 3, let w = d Xs ip. Then, by differentiating (jl.ip and (jl.2p with 
respect to x s , we have 

J2d Xi (A ij 9 Xj w) = 0, (4.1) 

where 

A ij = p(|V^| 2 )^- + 2p'(\V<p\ 2 )d Xi <pd Xj <p, (4.2) 

5u = 1, and 5ij = when i ^ j for i,j = 1, 2, 3. 

By (H2), this is a uniformly elliptic equation and A lJ are uniformly 
bounded (independent on k). 

Now, for large k > ko such that D k lies in the subsonic region, consider 
equation (14. lh in D' k . 

For any point P € rn9D£., by standard localized flattening and reflection 
arguments (here we use the Neumann boundary condition), in a ball-like 
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neighborhood B^ e (P) of P with the radius 3e depending only on 0, we can 
get by Theorem 8.24 in [12] that 

IIHI C «(B^(P)) < C \H\ L \B 3e (P)) < C '- ( 4 - 3 ) 

Note here that C, C, and a G (0,1) are independent of 93 and k. The second 
inequality follows from the fact that w is bounded. 

Now we analyze equation (|4.ip whose coefficients, after extension as above, 
satisfy 

||^|| C a (Bae (P))<ir, (4.4) 

for a constant K independent of k > ko- Therefore, by Theorem 8.32 in 
[12] . we see 

\\ w \\ch«(B e (P)) < C HHIc°(B 2e(P) ) < C "- ( 4 - 5 ) 

The constants C and C" are independent of k, and the second inequality 
follows also from the boundedness of |Vy|. 

Since T n dD' k is compact, the number J with V n C Uj =l B e (Pj) for 
ij G r n dD' k is independent of k. Then we obtain a uniform boundary 
Holder estimate of gradient of w with C independent of k: 

IMIci,<*(D fc n(u/ =1 B e (P.;)) ^ C - ( 4 - 6 ) 

Step 3. Uniformly global Holder estimate of gradient and second order 
derivatives of (p. 

Choose open sets D° k d D% g D fc such that D k - D° k C Uf =1 B e (Pj). The 

elliptic interior estimate (cf. Theorem 8.24 in [12]) to equation (|4.ip in D k 
tells us that 

\M\ ca .{p5 ) <c\\w\\ Loa p h) <c' 

with a' G (0, 1) and C, C > independent of /c. Without loss of generality, 
we assume that a' < a. 

Then we see that \\A i:> \\ ,,^77, < K holds for K independent of k. Now 

utilizing Theorem 8.32 in [T2] as above, we have 

lkll cl v ( ^o) < C \\w\\ Lx{Dk) < C" (4.7) 

with C" independent of k. 

Combining (|4.6|) with (|4.7|) . we conclude \\D 2 ip\\ Cct '^s < C with C > 
independent of k. This completes the proof of Theorem 1.3. 
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